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I.  FOREWORD. 

The  work  done  during  the  oeriod  covered  by  the  present  Final  Report 
constitutes  extensions  and  generalizations  of  the  results  found  under 
previous  grants  and  related  to  the  application  of  functional  analysis  to 
Fluid  ’’e  c  a  n  i  c  s  . 

■^he  specific  subject  dealt  viith  is  that  of  closed  splines  which  were 
introduced  ■f'or  the  first  time  by  the  principal  investigator  as  the  appro¬ 
priate  splines  to  solve  interpolation  problems  commected  with  the  flow 
fields  around  airfoils. 

This  work  has  lead  to  the  followino  paoers: 

1.  L.T,Napolitano;A  new  Characterization  of  Closed  Splines,  accepted 
for  publication  in  the  Italian  Journal :Aerotecnica  Missili  e  Spazio. 

2.  L.G.Napolitano:" Closed  Smoothing  Splines." 

The  second  oaoer  will  be  subtioned  either  to  the  previously  mentioned 
Italian  Journal  or  to  an  International  Journal. 

At  the  time  of  the  writing  of  the  present  Final  Report  the  final  choice 
had  not  been  made  yet. 

The  subject  Final  Report  contains,  in  extenso,  the  above  two  paoers 
which  constitute  Part  I  and  II,  respectively,  of  the  report  itself. 


PART  I 


A  new  Characterization  of  Closed  Splines. 


-3- 


A  NEV!  CHARACTERIZATION  OF  CLOSED  SPLINES 


I.  INTRODUCTION 

The  author  has  oreviously  defined  and  studied  new  classes  of 
spl’re  functions,  referred  to  as  closed  sd1 i nes  .  which  are  the  oroper 
splines  to  be  used  when  interpolating  data  prescribed  on  a  set  of 
points  belonoina  to  a  closed  contour  C  [ij  . 

The  theory  of  closed  solines  was  based  on  the  Hilbert-space 
aooroach  [pjand  was  thus  formulated  in  general  terms.  Different 
classes  of  closed  solines  can  be  obtained  by  approoriate  choices 
of  the  Hilbert  spaces  and  of  the  operators  acting  on  them. 

Two  anolications  of  the  general  theory  have  been  reported  in 
!_!  1  and  l_2  J  .  [ij  deals  with  closed  splines  interpolating  values 
of  a  function  orescribed  at  a  given  set  of  points,  [sj  deals  with 
Hermite  closed  solines,  i.e.  solines  interpolating  data  representing 
values  of  a  function  and  of  its  derivatives  prescribed  at  a  given 
set  of  points. 

Many  other  classes  of  closed  spline  functions  can  be  considered 
which  solve  other  interpolation  problems  of  practical  relevance. 

Thus ,  for  i nstance ; 


i)  The  orescribed  data  may  represent  a  linear  combination: 


ii)0ften,  exoecially  in  the  case  of  exoerimental  measurements, 
one  does  not  know  the  value  of  a  function  at  one  point  but. 


rather,  its  averaoe  value  over  a  small  intefiVal.  Thus  the  pre¬ 
scribed  data  may  be  of  the  form 


where  1  a i ,  bi)  are  the  subintervals  to  which  average  values  corre- 
SDond  and  the  weighting  factors  Hi  account  for  the  possible 
non - un i f ormi ty  of  measurements. 

iii) The  prescribed  data  may  represent  values  of  a  function  and  of 
any  number  of  its  subsequent  derivatives  or,  also,  a  linear 
combi  nation  of  them  . 

iv)  one  may  consider  trigonometric  spline  functions  (in  this  case 
what  is  being  is  not  previous  cases,  but  the  type  of  operator). 

This  li'St  is  only  indicative  and  is  noticeably  enlarged  if  two 
and  three  dimensional  solines  are  considered. 

For  all  cases  mentioned  above  the  characterization  and  proper¬ 
ties  of  classical  splines  (  i.e.  splines  defined  over  a  finite  in¬ 
terval  )are  well  established. 

Characterization  and  prooerties  of  closed  splines  (  i.e.  spli¬ 
nes  defined  over  a  closed  contour),  on  the  contrary,  are  established 
only  in  the  two  previously  mentioned  cases  of  normal  [l]  and  Hermite 
[3]  splines.  All  other  classes  of  splines  remain  to  be  characteri¬ 
zed  and  studied. 

The  aoolication  of  the  oeneral  theory,  as  formulated  in  [l ]  , 
entails  a  number  of  lenath^develooments  as  esemolified  by  the  contents 
of  [1 ]  and  [3]  .  It  appearred  thus  desirable  to  further  examine  the 
gene’^al  theory  with  the  aim  of  reformulating  it  so  as  to  make  it  pos¬ 
sible  to  utilize,  to  the  laeoest  possible  extent,  whatever  is  already 
known,  and  available,  from  the  corresoondi ng  classes  of  classical 
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s  D 1  i  n  e  s  . 

This  task  would  also  helo  to  shed  further  light  on  the  nature  and 
tyoes  of  differences  existing  between  classical  and  closed  splines. 

This  further  development  of  the  general  theory  of  closed  splines 
is  oresented  in  this  oaoer  for  the  case  in  which  the  space  X  is  an 
arbitrary  space  of  functions  defined  over  a  closed  interval. 

The  new  formulation  is  first  anolied  to  the  cases  already  studied 
in  [l  ]  ,  [  3  ]  i  n  order  to  further  clarify  the  essence  and  actual  pro¬ 
cedure  of  the  new  aporoach. 

Then,  as  a  futher  examole  of  aoolication,  the  new  class  of  closed 
spline  of  the  tyoe  (I)  above  is  characterized  and  its  properties  con¬ 
cisely  derived.  This  case  shows  the  extreme  usefulness  of  the  new 
aporoach  and  is  indicative  of  its  application  to  the  other  classes 
of  closed  splines  mentioned  above. 


2.  RELEVANT  RESULTS  FROM  ABS RACT-SPACE  THEORY. 

A  number  of  basic  results  of  the  Hi  1 bert-soace  spline  theory  are 
needed.  They  are  concisely  derived  in  this  section. 

Let  X,Y  and  Z  be  three  real  Hilbert  soaces;  T;  X  — > Y ;  A:X  — >Z 
two  linear  onto  continuous  operator  and  denote  by  <  , 
the  inner  product  and  the  orthogonal  complement  of  the  space  W ; 
and  by  R(B),  N(R),  s'the  range,  the  null  space  and  the  adjoint  of 
the  ooerator  B.  R(T)  and  R(A)  are  closed  in  Y  and  Z,  respectively. 

Definition  and  prooerties  of  the  spline  space  S  CX  correspon- 
dinn  to  T  and  A  hinne  on  the  following  egualities: 


i  e  S> 


=  <  .A,  /?<>_=  < 


X 


e  X 


X  e  Z 


(2.1) 
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where  t'  second  one  follows  from  the  definition  of  adjoint  operator. 
From  eqs.  (2.1)  it  follows  that; 

_  /  /■ —  \ 

(2.2) 


2)  I  x>^  =0 


X  G  hJCl  ) 

K  €.  NC(^) 


(2.3) 


Conseouentlv:  i_)  if  N(T1+  N{A)  is  closed  in  X  there  will  always  be 
at  least  one  X  G  ^  satisfying  eq.(2.I)  for  any  Ax  ^Z'(exi  stence 
theorem);  ii)  if  N(T)  /I  N(A)  =  |o]  this  /i  is  unique  for  any  given 
Ax  =  z6  Z.  (Uniqueness  theorem)  [zl  . 

The  equation: 


Ax=z  (2.4) 

with  z  an  arbitrary  but  fixed  element  of  Z  characterizes  the  "constra¬ 
ints 


Eq.  (2.3)  characterizes  the  spline  space  S  C  X  associated  with  the 
operators  ('T,A). 

Tt' 

When  the  space  Z  is  finite:  ,  with  its  usual  inner  product, 

the  operator  A  can  be  represented  as  a  set  of  (n)  linear  and  conti¬ 
nuous  operators  Ki  on  X: 


Ax= 


- y  Kv/  X  >^J 

Z  =  - 


(2.5) 


and : 


Yty  — 


<;  <  A,  C-fr 

A 


A .  <  IT.  X> 


R(A)=  Z  _  ^ 

R  (A  '  )=  /.N(A)]  =  a'  (Z) 


Furthermore  : 
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and ; 


/ 

dim  R(Ai=  dim  R(A  )=n 
co-dim  N(A)-  finite 


dim  q  =  finite 


(2.6) 


The  ^imteness  of  ‘J(T)  follows  from  the  uniqueness  condition  N{T) 
/n^N('A)=  -j^O^and  the  fact  that  N(A)  has  finite  co-dimension. 

From  the  well  known  relation  between  the  dimensions  of  subspaces 

i.  j_ 

dim!'.'  n'.'U  dim  W  .  dim  (  W,  f\  W,  )-dim  W, 

1  1  ±  i  -o 


and  tne  relations  holdino  for  a  linear  continuous  operator  with 
closed  rannes:  _  j_  j_  .  . 

H  fon„„s  -.Kat,  fo.  ft  : 

J:  n.  [H  fA'J  n  Rcr'j  ]  -  n  l^cr)]- 

-  Jl'ya.ncO  - 


(2.7) 


From  e  a  .  .  3  : 

.1  j'Ti  x>^^o  K  e  N(fi) 

hence  — Cl  RC*  /  so  that: 

T^e  ^  (r'J  n  R.  Cf^'J 

from  whicn;  uDOn  eq.(2.6) 

I  s:.  rt 

'  'U  'yyi. 


(2.8) 


REMARK  I 


In  essence,  n  is  the  number  of  constraints  (2.4)  which  are  imposed: 
i.e.  the  number  of  data  that  are  orescribed. 

REMARK2 

"he  dimension  of  T(S)  is  enual  to  the  number  (n)  of  data  prescribed 

i-n 

less  the  dimension  (q)  of  the  null  soace  of  / 

RE ‘-’ARK  R 

The  dimension  of  the  soline  snace  S  is  equal  to  the  number  of  data 
tnat  are  orescribed, 

REMARK  A 

All  above  is  quite  General,  '/'hether  one  considers  classical  or  closed 
splines  deoends  on  the  set  over  which  the  function  space  X  is  defined. 


3.  CHARACTERIZATION  OF  SPLINES. 

The  actual  procedure  most  often  used  to  characterize  classical 
snlines  of  a  niven  class  can  be  summarized  as  follows. 

Oiiven  the  In)  ooerators  Ki  one  finds  the  function  ^ 

sol vi nn  eauati on  (2.1) 


by 


<  'V<>, 


K' 


X 


(3.1) 


subject  to  the  conditions  of  eq.  (2.2): 


J 


X  e  ‘J 


(3.2) 


which,  in  neneral,  are  readi''  '  expressed  as  costraints  on  theAj. 
The  number  of  these  constraints  is  equal  to  dim  N(T),  i.e.  to  (q) 
ea.''2.6)J  .  Thus,  as  required,  by  eq.,(2.8),  ^contains  (n-q) 
arbitrary  paramenters.  constitute,  toqether  with  the  (q)  para 

meters  invo’ved  in  aettina  from  /  -4  to  4  ,  the  n  indipen- 

dent  parameters  upon  which  s  depends  ^  as  required  by  eq.(2.8)2 
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The  solution  of  eq.  (3.1)  can  often  be  put  in  the  form  /^2j  : 

t  ‘  .1 

-  c;3.2.; -  / Xe.  (3-2) 

where  the^Y  satisfy  the  (q)  constrants /  qljTs tions  to  be  addressed  are: 

Can  the  known  elements  ^enterina  the  characterization  of  classical 

'  c- 

splines  be  used  to  characterize  the  closed  splines?  if  so,  how? 

4.  CORRESPnNnENCE  BETWEEN  CLASSICAL  AND  CLOSED  SPLINES. 

The  problem  formulated  at  ^the  end  of  the  preceeding  section  is 
p r 0 oe r 1 CO r r e s 00 n de n ce  between  the  spaces  ana  operators  correspon- 
dina  to  the  classical  splines  and  those  corresponding  to  the  closed 
splines. 

In  this  section,  all  Quantities  pertaining  to  closed  splines  will 
be  denoted  with  a  subscript  (c). 

It  will  be  assumed  that  the  space  X  refers  to  functions  defined 
over  the  interval  [6,1  J  . 

Generalization  to  other  cases  may,  in  principle,  be  possible  but 
will  not  be  coside*"ed  here. 

The  set  of  points  at  which  data  are  prescribed  can  be  characterized 
by  the  values  Z^of  a  non-dimensional  pa rameterZ" .  These  values  are 
such  that ; 

(9  =  r,  <  <t - <  ^ 

For  closed  splines,  the  ' s  are  the  curvilinear  coordinates  along 
an  arbitrary  contour^  C,  sufficiently  smooth  and  regular,  measured 
from  one  of  its  points  (  f,  -0)  and  normalized  with  respect  to  its 
length  =I).The  common  point  of  '‘coordinates  ^='^^^Xwill  be 

referred  t  as  closure  point. 
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A  first  corresDondence  is  thus  established,  for  both  X  and  X 

c 

are  defined,  unless  an  inessential  stretchino,  over  the  same  set 

IJ  . 

The  elements  of  must  satisfy  continuity  requirements 

at  the  closure  ooint  whose  number  and  nature  depends  on  the  type 
of  soace  X.  Hence; 

X  C  X 
c 

The  last  requirement  for  the  oroblem  at  hand  to  be  well  posed  is 
that  the  classes  of  splines,  be  the  same.  This  means  that  ^  ^ 

^  Z  and  the  ooerators  Tc,  Ac;  Tc:  Xc— i»Y;  Ac;  Xc-^Z, 

be  the  restrictions  of  the  operators  Z  and  A  to  Xc. 

The  most  relevant  consequences  of  these  correspondeces  will  now 
be  anal yzed . 

To  beoin  with;  from  Xc  CX  it  may  follow  that  N(Tc)  C  N(T)  and: 
d.i'VVU  urn 

so  that,  if  Sc C Xc  C  X  is  the  space  of  closed  spline  functions; 

‘12(SJ=  ft  w - TCsj 

Furthermore,  when  there  are  also  differences  in  "nature" 

between  and  ^  oecause  they  satisfy  different  constraints. 

To  elucidate  this  point  consider  the  set  as  a  base  and 

denote  by  F C  Y  the  subsoace  spaimed  by  this  base.  On  the  olansi- 
ble  assumption  that  the  sAd  's  are  indipendent,  F  is  an  (n  +  I)- 
dimensional  soace.  Consider  now  the  element  f  <£*  F  defined  by: 


T  T 

-  1  1  - 


('  ^ 

where  thi='  -^are  now  arbitrary.  Clearly  TSC'F  and,  wore  specifi¬ 
ed  11  ^  ^  the  satisfy  the  (q)  relations: 

Z  f.  <  K-  X  >  -  0  X  e  U(^r) 

/''‘i  y  »✓ 


If  N(ic)CN(y)  and  the  satisfy  only  the  ^  ^  relations 

Yi-rl  -  _ 

Z  A  -  <  x>^  X  eMCu) 

i..  -  X" 


(4.2) 


f  beionns  neither  toT(Sjnor  to  Tc(ScJ 

The  difference  between  ^and  is  now  evident.  It  is  solely  due 

‘c.  '  - 

to  the  smaller^of  constraints  imoosed  on  the^^jby  the  smaller  dimen¬ 
sion  of  M(Tc)  and,  does  not  involve  the  Hence  it  must  refledt 

the  different  behaviour  ofy^andy^  at  the  closure  point.  Foryy  the 
closure  Doint  can  be  any  noint  and  thus  the  features  of  ^  at  the  clo¬ 
sure  Doint  are  the  same  as  those  at  any  other  point.  It  will  be  pre¬ 
sently  seen  that  this  criterion  leads  to  the  new  characterization  of 
yy  Which  makes  uses  of  all  known  orooerties  of  ^  and 

Embed  a  subsoace  od  F  into  an  higher  dimensional  space  Fc  such  that 
Tc(_ScJCFc. 

Only  p1  ptnents  </^or  ^  of  the  basis  of  F  is  rele- 

vant,  as  it  can  be  inferred  from  the  continuity  conditions  at  the  clo¬ 


sure  point. 

Thus  the  base  of  Fc  contains  the  first  n  element  .  On  the  other  hand, 
in  order  for  Tc(^ScJd;Fc,  the  dimension  o  of  Fc  must  be  such  that  once 
the  qc  constraints  (4.2)  and  the  remaining  (  q-qc)  satisfied  by  ^ 
but  not  lU  ?ire  endorsed,  one  obtains  a  (n  +  I-qc)  dimensional  space,  as 
requ i red  by  eq .  (4.1). 

Hence  p- q c- ( q -q c ) =n +  I -qc  ,  i.e.  p  =  n  +  ( q-qc- 1 ) . 


-I?- 


The  element  ^  ^  ^  can  thus  be  formally  written  as: 

-  « 

I  -  2  A..  ^  f-  f' 


(4.3) 


where  the  must  be  such  that  TcScCFc  and  can  usually  be  found, 
‘d 

by  insoection,  once  the  expression  of  the  ^.iis  known. 

The  new  formulation  for  the  characterization  of  closed  spline 
is  obtained  by  combininn  eq.(4.3)  with  the  oreviously  mentioned 
c  r  i  t  e  r  i  0  n  . 

Let  S  any  classical  spline  space  defined  over  the  interval 
L0,lJ  and  call  interior  point  Pi  any  point  within  this  interval 
and  extremal  the  left  and  right  extreme  of  the  interval. 

The  characterization  of  an  element  s  fc' S  defines  the  func¬ 
tion  y^and  the  properties  of  Ts  at  interior  points.  In  particu- 
1 ar  one  knows : 


I)  the  number  and  types  of  discontinuities  at  interior  points. 


2)  the  number  and  types  of  continuity  properties  at  the  interior 
DO i n ts  . 

Denote  by  the  set  of  interior  point  discontinuities  and 

by  |C  j  the  set  of  interior  point  continuity  prooerties. 

Then  the  following  characterization  of  closed  splines  of  the 
same  class  holds: 

New  characterization  of  closed  splines. 

For  any  class  of  classical  solines  defined  over  a  closed  interval 
the  corresnondi no  closed  spline  can  be  expressed  as: 
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where  is  niven  by  the  classical  soline. 

inspection  and  the  (n+q-q  +  I  )  coefficients 
minedj^imposino:  at  the  closure  ooint: 


is  inferred  by 
^'and  ^  •  are  deter- 


I)  the  connruence  of  the  di  sconti  nui  tes  of  the  set  ^  ■ 

Z]  the  vanishino  of  the  discontinuities  corresponding  to  the  set 

■ 

Condition  (I)  can  be  expressed  for  an  arbitrary  function  h  discon 
tinuous  at  the  points  Pi: 
as  (  see  /.  3  7  for  the  proof): 

Z  /  4-  i  (4. 

i 

where: 

K:k- 

remark  I 

Oncey^  pc  sc  is  known,  obtaining  sc  is  a  trivial  matter. 

'c 

REMARK  2 

For  classical  spl i nes^condi tions  (2)  are,  by  definition,  not 


satisfied. 
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5.  EXAMPLES  OF  APPLICATION  OF  THE  NEW  FORMULATION. 

The  classes  of  splines  for  which  the  closed  ones  have  already 
been  derived  /_  I  ]  ,  /^ojwill  be  dealt  with  first  to  help  clarifying 
the  procedure  to  be  used. 


5  .  I  Normal  sol  i  nes  [  I  J  . 

Soaces  and  ooerators  are  chosen  as; 


X=  ft 


o 


d-h  1 


r=  Al-  =  C). 


A 


a 

:  6  [i, 


where  is  the  Hilbert  soace  of  real  functions  having  sjquare- 

interarable  derivatives. 

I 

From  [2  ]  ouQ  gets 


i) 


f- 


.  j  L  £  n. 


(3- 


I  i 


with: 


c 


-  2: 


-'i  - 


(z-z:) 
0 


if  (z-rc.)fo 


if  {v-r:)<LO 


1  1 


a-’-  -  ^ 


^  I  n 


iii)the  derivative  of  order  (q-I)  of  i%  discontinuous  (set  {D^ 
of  dimension  I) 


iv)  and  its  derivatives  up  to  order  (q-2)  are  continuous  (set 
•(Cjof  dimension  q-I) 
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For  closed  splines  defined  over  the  closed  contour  C  dim  N(Tc)=dim 
NH^{C)  =  I  [x],henceq^=I. 

Thus,  upon  eq.(4.4) 


0 


o 


with  the  (r+q)  coefficients  satisfying  the  conditions: 

IL. 


Z  i: 


i.  -  / 


(JlJ 

-r  ci')-o 


(5.1) 

i:  fee?  ^ 


Remarks 

For  conventional  splines: 

.  The  conditions  analogous  to  the  first  of  eqs.(5.I) 


Tt 


(5.2) 


expresses  the  discontinuity  at  the  point  Pn  +  I  in  terms  of 

those  {  A  ■  )  at  all  other  points. 

The  set  of  equations  due  to  the  constraint  x6  uC  leads  to 
the  conditions: 


(5.i) 


.  Equation  (5.2)  is  also  a  congruence  condition  since  outside  the 
i nterva 1  [3,1  J  is  zero.  Indeed; 
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from  which  eq.(5.I)  follows  since  j  is  constant  in  each  sub¬ 
interval:  J 

5  .  2  Herrrute  splines  ^3  "] 


Soace  and  ooerators  are  now  chosen  as: 

x-if'  ;  '-.1-//"  . 

A  A 

so  that  now  dimZ=^«^(see  remark  I, Sect. 2). 

From  lone  nets:  ^  n  >*  7 

„  f,.jt-r,jr/crof ;  ^‘-Cr-=A /ft-t).' 

iii)  The  derivatives  of  order  (q-I)  and  (q-2)  are  discontinuous 
(set  -[oj  of  dimension  2). 

iv) y^and  its  derivatives  of  order  up  to  (q-3)  are  continuous  (set 

of  dimension  (q-2)  . 

m 

Thus  ,  upon  eq  .  ( 4 . 4) : 


•A  /t  u  ^  X 

K 


r  I 

Ut  ■  A/ 


,  I'iy 

T  =  X  ^ 


ir  g  si  » 2-  /  _  J  .  J 
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The  (2n+q)  coefficients  satisfy  the  q  equations: 

'K-  - 


I 


^  A  '  ~o  y  R  f 

c^,  ^  '^9-/ 


The  second  equation  follows  f  ron’^  eq  .  (  4 . 5  )  [s]  . 

For  conventional  splines  the  satisfy  the  (q)  conditions  fz'] 


nil  __ 

z 


C-  I 


which  lead  to  the  same  conditions  (5.3)  at  the  extremal  points. 
The  space  of  Hermite  closed  splines  is  of  di mens i ons ( 2n )  . 

One  more  class  of  closed  splines  will  now  be  considered.lt  is 
new  insofar  as  it  has  not  been  introduced  before.  Its  classical 
counterparts  is  well  known.  The  new  formulation  will  be  used. 
Results  are  stated  concisely. 


5 . 3  Data  are  linear  combinations  of  values  of  a  function  and  its 
first  derivative. 

The  spaces  and  operators  are  defined  as: 
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whsre  the  constants  vanishing.  From^sJ  : 

i) 


.  o',  ^Z'-^O-r  _  p.  (^- 

ii)  As  the  operator  7^  is  the  same  as  in  the  previous  two  cases 


the  be  the  same. 

7 

iii) The  derivatives  of  order  (q-I)  and  (q-2)  are  discontinuous 
(set  |oj  of  dimension  2)  and  are  related  by  the: 

iv)  y' and  its  derivatives  up  to  order  (q-3)  are  continuous  (set 
.j^cjof  dimension  (q-2)). 

Thus  ,  upon  eo  .  (  4 . 4  ) 

where  the  (n+q)  the  coefficientsR,J,'sat1sfy  the  equations: 

V  ‘ 

'I:  [o^c  z'i 


^  I,-*-" 


c  -  • 


{^] 


c  ^'"7 


O 


j£  k  ^ 


The  relevant  oroperties  of  these  classes  of  closed  splines  are: 
I)  ^is  a  Dolynomital  of  degree  (2a-I)  in  each  open  interval 
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2)  i  s  continuous  on  C  toaethen  with  its  first  (2q-3)  derivatives 

3)  at  each  ooint  ( Zj )  the  discontinuities  of  the  (2q-2)-th  and  (2q-I) 
-th  derivatives  are  such  that: 


S 


s:.  + 


y  L  6  U, 


4)  ^is  such  that: 

The  set  of  functions  of  //  (C)  satisfying  the  first  three  conditions 
contitutes  the  n-dimensional  closed  spline  space  Sc. 

The  uniqueness  conditions  read  [Sj  : 


^^Cz)-o  "  . 


y  I  e  [', «  J 
C 


imoly  k(Z  )  =  O 


6 .  Concluding  remarks . 

A  new  formulation  for  the  characterization  of  closed  splines  has 
been  derived. 

This  new  formulation  makes  it  possible  to  use  a  number  of  facts 
already  known  from  the  theory  of  classical  splines  and  helps  gaining 
a  deeper  understanding  of  the  differences  between  the  two  types  of 
splines. 

The  aoolicatio  of  this  new  f ormu lation  is  esemplified  for  three 
classes  of  splines.  Two  of  them  had  already  been  derived  J, 
/_2]with  the  direct  approach.  They  hav’e  been  considered  in  order  to 
further  clarify  the  essence  of  the  new  formulation. 

The  third  class  has  been  deduced  here  for  the  first  time. 
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The  new  formulation  by  passes  lono  and  tedious  devel opments and 
makes  it  almost  immediate  to  deduce  the  closed  spline  counterpart 
of  known  classical  splines. 

The  new  formulation  concerns  only  the  characterization  of  closed 
splines.  All  other  relevant  properties  (  exi stence^ uniqueness ; 
extremal  properties  eouivalent^with  a  minimum  problem)  are  those 
expounded  and  discussed  in  the  original  derivation  of  the  theory 
I  3  and  have  consequently  not  been  repeated  here. 
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CLOSED  SMOOTHHlTi  SPLINES 

I.  INTRODUTION 


In  comoutationsl  aerodynamicsit  is  often  necessary  to  solve  Interpol 
tino  oroblems  related  to  airfoils,  i.c.  to  closed  curves. 

Usino  classical  spline  functions  to  interpolate  airfoil's  ordir.ates 
prescribed  on  a  finite  set  of  points  is  unsatisfactory  on  many 
imoortant  acconuts,  as  fully  discussed  in  [l]  .  Similar  situations 
arisei  in  all  other  interpolation  problems.  The  case  of  Hermite 
£ni'’'ne  functions  (  the  data  to  be  interpolated  represent  va''ues  of 
the  ^'unctions  and  of  its  first  derivate  at  a  given  set  of  points) 
was  discussed  in  C2 ]  • 

The  main  short  cominos  can  be  briefly  described,  in  general  terms, 
as  follows,  '.'hen  the  data  to  be  interpolated  are  prescribed  on  a  set  of 
ocintsevjbelonaing  to  a  closed  curve  (e.g.  the  airfoil)  and 
conventional  snlines  are  used,  then  the  first  and  last  points  must  be 
taken  as  coincident  (closure  point-usually  the  airfoil's 

trail  inn  edce).  As  a  consequence: 

ilunwanted  (and  unnecessary)  discontinuities  are  introduced  at  the 
closure  noint; 

ii'ia  number  of  derivaties  of  the  unterpolatino  spline  vanish  at  the 
closure  noint  and  this  may  turn  out  to  be  too  penalizing  *n  practical 
a  on  1 i ca ti ons ,  especially  with  splines  of  low  order; 
iii)the  decree  of  smoothness  of  the  interpolating  function  is  neither 
ouantifi cable  nor  uniform  over  the  closed  contur  with  the  accuracy 
be  inn  very  poor  at  the  closure  point; 

iiii'the  interpolation  function  does  not  satisfy  any  minimization 
problem  and  thus  all  the  ad  vantages 'comec ted  with  this  fact  are 
lost. 

To  overcome  all  these  shortcomings  the  author  developed  a  theory 
for  a  new  classes  of  solines,  referred  to  as  closed  splines  cu. 
based  on  the  general  abstract-space  spline  theory  detailed  in  fS} . 
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"he  aost'^act  soace  aoDroach  presents  the  following  main  advantages: 
i;  existence  and  uniqueness  can  be  proved;  ii)general  characteri= 
zations  of  the  splines  can  be  given;  iii)  the  equivalence  between 
an  i n te r 00 1  a t i on  problem  and  a  minimization  problem  can  be 
estaolisned;  iiii)  extremal  and  other  important  properties  of  the 
soiines  can  be  demo  seated. 

This  oaoer  addresses,  with  the  same  approch,  another  important  class 
of  problems  of  frequent  occurrence  in  computational  aerodynamic^: 
namely  the  one  solved  by  smooth! no  soli nes . 

To  clarify  the  issues  involved,  consider  the  case  in  which  the 
data  to  be  interoolated  represent  the  values  of  a  function  at 
a  niven  set  of  points  (Qj. 

In  many  cases  the  data  prescribed  are  not  exact  but  approximate, 
i.e.  affected  by  either  errors  or  uncertainties.  It  thus  makes 
little  sense  to  have  the  interpolating  function  f  assume  exactly  the 
values  V  3t  the  points  and  one  would  rather  like  to  achieve  a 

ii-L 

seEtable  compromise  between  the  "smoothmess"  of  the  interpolating 
function  and  the  "approximation"  of  the  prescribed  data 
The  classes  of  splines  achieving  this  compromise  are  called 
" 5  moo  t  h i n  n "  splines. 

The  formulation  uf  smoothing  problems  solved  by  nlassical 
imoothinq  splines  and  their  properties  are  well  known  ["3^  and  can  be 
summarized  as  follows. 

Given  a  closed  interval,  reduced  through  stM^able  normalization 
to  /.  iV;  fj;  (wr/)  points  (2^(  rj)  with: 

^  r;  c  -  ^  ^  ^ 


and  real  numbers/^'(  the  smoothing  spline  function 

of  order  'j  £  (  n  -t-  i  ]  cirresp  ending  to  the  set|'2^^^^^],  where^is  a 
positive  constant,  is  the  unique  function  0“£ /f  ^7  wh  i  ch  solves  the 
following  minimum  problem: 


Jm/yu 

P.- 


{ j'[ 


(I 


-24- 


where  I's  S  the  Hilbert  space  of  real  functions  ^[z)  defined  on 
[^]i]  and  having  a  square-intergrabl  e  ^/i(deri  vati  ve  The  space  // 
quantifies  the  degree  of  smoothness,  the  coeff i ci ent  y characteri zes 
the  relative  importance  that  one  assignes  to  the  smoothness  and  the 
a Doro X  i  ma t  i  on  of  data. (More  generally,  different  weights  y^can  be 
prescribed  for  each  point). 

When  one  sets  f(2(;)  =i(;the  minimum  problem  (I. I)  reduces  to  the 
one  pertaining  to  the  interpolating  splines  il,3j  and,  for  q=(«+/  ) 
the  interoolating  spline  is  the  unique  poli  nominal  of  degree  {  n  -t  !  ) 
passing  througlithe  points 

Smoothing  and  interpolating  splines  belong  to  the  same  subspace 
S  <- //  of  real  functions  s  (tO  defined  on  and  such  that  ['ij: 

a)  s('r)  is  a  polynomial  of  degree  (2q-I)  in  each  of  the  open  intervals 

r:  fij  [1  ^  ^  ^ 


b)  s(  Z  )  and  its  first  (2q-2)  derivatives  are  continu  e«-s  o  n  i  7 

c)  the  derivatives  of  s{Z)  from  orders  q  to  (2q-I),  included, 

vanish  at  the  points  and  (j^i  '  0- 

Ti)e  interpolating  spline  for  the  set  the  unique  element  (T 

of  S  such  that;  O'  j  =  j  i  t  i' , 

The  smoothing  spline  (Tfor  the  set  the  unique  element  of  S 

ft:  -  ^Lzc)]  - 


where  : 


is  the  discontinuity  (jump)  of  the  (2q-I)-th  derivative  ofcrat  the  poin 

Q.  I  - 

Since  the  smoothing  spline  belongs  to  the  same  space  S  as  the 
interpolating  spline,  it  is  evident  that  when  the  set  of  points  belong 
to  a  closed  contour,  smoothing  problems  exhibeit  all  the  previously 
mentioned  shortcomings. 

Thus  the  need  for  developing  a  theory  of  closed  smoowOlng 
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splines  exists  and  is  substantiated  by  the  same  consideration 
that  were  made  in/j/^zjwhen  dealing  with  interpolation  problems 
on  closed  curves . 

In  this  paper  the  general  theory  of  closed  smoothing  splines  is 
developed.  It  forms  the  basis  for  the  subsequent  consruction  and 
study  of  different  classes  of  smoothing  splines(  normal,  Hermite, 
Fou  r  i  e  r  ,  Tr  i  gonome  t  ri  c  ,  local  averange  and  so  forthTzJ). 

As  a  first  example  of  application  of  the  general  theory,  the 
normal  closed  smoothing  splines,  corri spondi ng  to  the  problem 
discussed  above,  are  derived  and  their  propertes  studied. 

As  in  Li]  and  ,  the  presentation  does  not  follow  the  logical 
order.  To  facilitate  a  more  w^e  spread  comprehension  and  use 
of  normal  closed,  smoothing  splines  their  definition, 
characterization  and  properties  are  first  stated  without  proof 
in  section  (2),  The  general  Hilbert  space  theory  of  closed 
smoothing  splines  is  then  developed  in  sections  3. 

The  last  section  4  presents  the  proof  of  the  statements  made 
in  Section  2  concerning  existence,  uniqueness^ characterization, 
extremal  and  other  relevant  properties. 
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2. NORMAL  CLOSED  SMOOTHING  SPLINES. 

2.1  DEFINITION 

Let  C  be  a  sufficiently  smooth  and  regular  closed  contour. 
Denote  by^the  curvilinear  coordinate  along  C  measured  from  an 
arbitrary  initial  point  and  normalized  with  respect  to  the 
lengthof  C  so  that  0-'^-^.  The  pointfiwill  be  referred  to  as 
the  cl osure  poi nt  of  the  contour  and  it  is  characterized  by 
either  values  2r=.(9^and  I  of  . 

Consider  (n)  arbitrary  success’ <e  pointsQ*(  'K.  )  on 

C,  let  be  their  curvilinear  coordinates  with; 

D  -z,  -  <  i 

and  prescribe  n  real  numbers  (fej)  and  a  positive  constant 

The  closed  smoothing  spl  i ne  (r(^Jof  degree  (q)  corresponding 
to  the  (  n)  conples  {  Zu  )  is  defined  as  the 

unique  element  of//'(  C)  such  that: 


_ 


'ie/fVcJ 


r )]  ?  •■r;  /  /  ^ 


(2.1) 


In  order  for  <T'6ff(C)  to  be  such  a  closed  smoothing  spline 
it  is  necessary  and  sufficient  that: 

a) G'(%)be  a  oolinomial  of  degree  (2q-I)  in  each  open  interval: 

Co  -  z-^)  ^  — j  (zc  ,  ^  ^  d  ) 

b)  CTi^jbe  continuous  on  C  together  with  its  f^rst  (2q-2)  derivatives,  i.e 


y  k£r 


(2.2) 
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)  (5iz)  be  such  that; 

/  t;,-<rirc)J  <  0  '-> 


n 


where:  . 

^'^^r.ctes  the  uM' scCMti  r.u’ *-y  of  t^e  (2q-n-th  derivative  at  the 
point  Q.*,  the  di  scouti  nui  ty  at  the  closure  point  being  defined 


as  : 


The  set  of  functions  4  6/f(C)  satisfying  the  first  two  conditions 
constitute  the  space  -S^of  the  closed  spline  functions  correspon= 
ding  to  the  set  a  subspace  of  jfiC)  of  dimension  (n). 

Interpolating  and  smoothing  spline  functions  relative  to  the  set 
^ZjJare  subspace  of  the  same  space  ^  j 
fot  Interpolating  splines  aafee  reolace-^i^  the  conditions 

2.2  CHARACTERIZATION 

Introduce  the  function^sl: 

Z^~l 

^  if  (z-n)zo  ,2  4, 

(£"Wf  "  I  O  if  Z'-'Jr'  f  ^ 

An  element  56//(C)  belongs  to  the  subspace.^of  closed 
spline  functions  corresponding  to  {2ijif^it  is  representable 
as:  :  -K, 

i  (2-5) 

where  the  n  coe^f  i  ci  ent /4.and  theZq  coef  f  i  ci  ents  0  s ati  sfy  the 
following  (2q)  equations  /l7  ^  ^ 

o  '  k  ^  -I 


o 
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The  function  defined  by  eqs.  (2.4)  is  continuous  on  C  together 
with  its  first  (2q-2)  derivatives.  Its  (2q-I)-th  derivative  is 
discontinuous  at  the  discontinuity  being  equal  to  [3  ]  . 

The  (2q-I)-th  derivative  of  3 (j^)  is  given  by; 

p.,-' 

Hence  it  is  piecewise  constant  in  each  open  sub-interval  and: 

i 

X 

-V'-j- 

n)  coef f i ci entsX'represent  the  values  of  the 


1^’"'-  p. 


V 


1  n 


(} 


-/ 


Thus  the 


discontinuity  of  at  the  points  ifi: 

Since  for  it  is  2)- =  and  Z.  ~  1 

■P'e  closure  pointiP^is  thus  in  no-way  differentiated  from  the 
other  (interior)  points. 

Given  (n)  arbitrary  real  numbers  and  a  real  constant 

^  there  is  a  unique  closed  smoothing  spline  corresponding  to 
the  sets  f  such  that: 


y  fir;  - 


/i  ^ 


(2.8) 


Consequently  the  system  of  (n-t^)  equa  t  i  ons  (  2 . 6 )  ,  (  2 . 8)  , 
linear  inp,andy|*,  admits  of  a  unique  solution.  When  the  (  )  coef 

ficients  apoearing  the  eq  (2.5)  are  so  ob»ta  i  ned ,  becomes 

i  p-£(>udl  O'4'tJ  -vo-iUZ/TV^ 

the  unique  smooth i ngl^obl em  (Z.I). 

When  eqs.  (2’,S}  are  reolaced  by  the  equations: 

f£  c  ^  n. 

the  solution  of  the  system  formed  by  these  equations  and  eqs. (2.6) 

m 

is  still  unique. The  new  coefficients  ,  thus  derived,  when 
sostitu^d  in  eq .  (2.5),  define  the  'inique  interpolating  closed 

s  p  1  i  ne  ^(4jwh  i  ch  solves  the  minimum  problem[lJ'. 

•  ,cis->u 

formally  obtained  from  the  problem  (2.1)  by  setting  ^ 
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2.3.  PROPENTIES 

The  closed  smoothing  spline  iT  corresponding  to  the  set 
has  the  properties  a),b),c),  described  in  section(2.I). 

Plus,,  obviously,  all  properties  of  closed  splines.  Thus,  for 
instance,  the  first  of  equations  (2.6)  represents  the  congruence 

C  -r 

conditions  for  the  discontinuities  0^^'"  for  any  element  6 

and  any  function  ^6  H  (pj  the  following  relation  holds  £lj. 


L=-  / 


(2.9) 


In  parti  cul  ar ;  .  9/  ■,  / 

^  v^elo-l  feffVcJ/ 

I  I  L  I 

When  ^  is  a  closed  smoothing  spline  G^eq.  (2.3)  holds  and 
eq.  (2.9)  becomes:  ij-  - 

j  p  Jz  -  f  A  z 

and  eq.  (2.10)  still  holds  fors=  (T". 

The  minimum  problem  solved  by  G'(^)  is  a  particular  case  of  the 
following  more  general  extremal  properties. 

Given  ^n)  arbitrary  but  fixed  (7^),  if  6"  is  the  unique  element  of 
^such  that: 


y  f  =  C-0  ^ 


then  : 

a  )  for  any  f  <£  (  C )  : 


-0 


(2. II) 


/"c 

and  (T  i s  the  unique  element  having  this  propety. 

The  original  minimum  probl6»)t(2. 1 )  follows  from  eq.(2.I2) 


for  s=o. 
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2.4.  COMPARISON  WITH  CLASSICAL  SMOOTHING  SPLINES. 

The  classical  smoothing  spline  over  (n-t)  )  points  Qi^s  repre  = 

rentable  as  [3  I  ;  ,  -v  Zsi-f 

SiZ]--2i  - p—  (2.13) 

(^/  V  t*  (J.ij  -OJ 

The  {  /It  coefficients  represent  the  unique  solution  of  the 


followinc  system  of  equations: 

£  X  :o^k^~Ci-0 

.of  classical  splines  is  the  set  of  all  ^of  the  form 
given  by  eq.  (2.13)  with  the  coefficients  satisfying  only  the  first 
(q  )  eqs  .  (2 . 14) . 

When  eqs.  (2.14)  are  used  for  the  set  of  points  belonging  to  ^ 
closed  contour ,  (9,5  and  .  Since  the  coefficients^',// 

*ntl  f 

in  eq.  (2.13)  are  uniquely  determined  it  follows  that: 

i)  at  the  closure  point  (^.-^.^/  )  the  condition  guaraktees 

only  the  continuity  of  0~  (^J  :  its  derivatives,  up  to  the  order 
(^- 1 J i nc 1 uded  ,  will  not  be  continuous (uni  ess  the  problem’s  data 
satisfy  very  paricular  synmetry  conditions)  and  their  disconti= 
nuities  are  uniquely  determined. 

ii)  the  subsequent  derivatives,  up  to  the  order  (2q-2),  are 
continuous  but  vanish  identically  (see  sect.  I.) 

This  has  to  be  contrasted  with  the  closed  smoothing  spline 
which  is  continuous  throught  together  with  its  first  (2q-2) 
derivatives,  the  properties  at  the  closure  point  being  in  no-*^®-^ 
different  from  those  at  the  other  points. 

The  classical  smoothing  spline  has  (q4M-»0  parameters:  if  one 


introduces  (  ^  -/  )  additional  parameters  to  impose  the  missing 
continuity  of  the  )  derivatives  at  the  closure  point,  one 

would  get  a  function  with  a  total  )-r(^~/  )*(2q^»t') 

parameters.  These  are  indeed  the  number  of  parameters  appearing 
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in  the  definition,  eq.(2; 5),  of  the  closed  smoothing  spline.  This 
rather  naive  approach  does  not  however  indicates  how  the  additional 
oa rameters shoul d  be  introduced,  what  systemof  equations  should  be 
formulated  for  their  unique  determi na ti on ( cl  earl y  the  greatest  part 
of  eqs.  (2.14)  cannot  be  retained)  and,  most  important  of  all,  no 
hints  cnuld  be  obtained  as  to  how  to  determine  the  properties  of  these  n 
functi ons . 

The  obstract  approach,  leading  to  the  results  presented  in 
sections  2.1  through  2.3,  provides  in  the  more  natural  and  correct 
maimer  the  answers  to  all  questions  above. 
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3.  the  smoothing  spline  theory  in  abstract  space 


3j  arbitrary  helbert  spaces. 

0 

The  needed  basic  results  of  the  Hilbert  space  formulation  of 
smoothing  spline  function  theory  are  summarized  here  for  ready  and 
convenient  reference  (  see  [,3 1  for  greater  details). 

LetX,Y,Z,  be  three  real  Hi^lbert  spaces  with  norms  {{•  ii  ^ 

II'  'j  11  '  ;  and  let  I  ■  X^->  ;  /}-'  X  ~>  ^  be  two  1  inear 

continuous  ooerators  which,  without  any  loss  of  generality,  are 
supDosed  to  be  onto. 

Given  any  fixed  and  any  real  constant  <^>0  consider  the 

followino  minimum  problem  .  ;  ,,  ^  i,  7  \ 

II  T^/l^  tfl 

The  e  1  emen  t  i  f  it  exists,  is  called  the  smoothing  spline 
corresponding  to  (  7//^  5* )  • 

Introduce  the  Helbert  space  P,  cartesian  product  of  Y  and  Z 
endowed  with  the  scalar  product: 


(3,1) 


where  <,  >p  denotes  the  inner  product  in  the  Helbert  space  F. 

Then  if  i:  y->?is  the  linear  continues  operator  defined  by: 

--  CTy,  A  .(I  €  P  Ci-iJ 

and  if:  r  t  O 

10,27  e? 

where  0  denotes  the  null  element  of  Y,  the  minimum 


problem  (3W)  can  also  be  formulated  in  the  space  P  as: 

II  ^  ^  -  fo  //p 


a; 


where 

(3-2)  • 


denotes  the  norm  of  P  associated  with  scalar  product 
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The  followinn  existence  and  uniqueness  theorem  holds. 

THEOREM  I.  The  solution  of  problem  (3,/)  exists  for  any  z  ^ 
iff  in  closed  in  X  and  is  nuique  iff,  in  addition, 

M  CT  j  0  ~{^j.  Here  N(B)  denotes  the  null  space  of  the 

operator  B  . 

The  characterization  of  the  smoothing  spline  function  is  given 
in  the  following  theorem. 

THEOREM  2.  Under  the  hypothesis  of  theorem  I,  the 

smothing  spline  corresponding  to{T,A,Z,  9)  ^ 

<  la  ~  1=^^^  Ly  >2  -O  ^  x:  €:  X  (^3. 

Thus  the  smoothing  spline  space  S  is  the  subspace  of  X  defined  by: 


or 

L; 


equivalently,  on  account  of  the  definition  of  the  adjoint  of 
by: 


S  =1  deX/  j 


where  the  prime  denotes,  here  and  in  what  follows,  the  adjoint  of 
an  operator;  RCb*)  is  the  range  of  the  operator  B  and  the  adjoint 
L‘  of  L  is  defined  by: 

The  second  definition  follows  from  the  first  one  on  account  of 
the  oroperties  of  adjoint  operators  and  of  the  fact  that 

f>  R  ^  6  'R.  'J-' 

^  la-  <  I 

From  the  definition  (3,?)  it  also  follows  that: 


(3.9  ) 
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REMARK  I. 

It  Can  be  rea 1 y  shown 
thus  coincides  with  the 
the  operators  T  and  A. 


that  S  does  not  depend  on  )  and 
sqace^of  the  spline  fonctions  defined  by 


i  t 


THEOREM  3.  For  any  fixed  there  exists  a  ncique  element  G"6 

such  that : 


Pi^  t  ±  Qrr  ^  e 


—  r  "7 

):hi«'  space  S  can  be  characterized  as/_3j  :  a 

where  N  (  •  )  denotes  the  null  space.  Hence  ^  S  iff; 

r'--3  ffiCfl'i 


Hence: 


and,  upon  eq,  ( Sis’)  ,  S*'=  ^  since  ^ 

^  fi-'  ft  ^  ^ 


and 


denotes  the  inverse  of  an  o;  rator. 


REMARK  2. 

The  corresponding  theorem  for  interpolating  splines  ensures  the 
existence  of  a  unique  element  O” S>  such  that 

The  characterization  of  a  smoothing  spline  which  leads  to  its 
pratical  evaluation  in  contained  in  the  following  corollary  of 
Thorem  2. 

COROLLAPY  (  Character!  zati  onj  cff~J^oothi  ng  spline  corresponding  to 
(T,A,z,  j)  iff  there  exists  A  ^ sjjch  that: 

'  I  "3.  i 


Indeed,  from  eqs.  (3.2)  and  (3.5)* 
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L'Ls 

==  7'7'cr  =  ?  R'lz-fis-) 

from  wnich  eqs.  (.3./<)  follow  with  A=  ^  (z-O^)  . 

The  first  of  eqs . 

(3[/)  shows  that  A  €  t^(j)  ^  where  (  denotes  the 

orthoaonal  subsoace^and  the  second  equality  follows  from  the 
properties  of  adjoint  operators. 

The  estremal  properties  of  the  smoothing  splines  are  condensed 
in  the  following  theorem. 


THEOREM  If  z  is  an  arbytrary  but  fixed  element  of  .Z^  and  0“i  s  the 
unioue  elemente  of  S  satisfying  the  condition  (3i.(o)  ,  then: 
a  )  for  any 

I)  7  (,5"-  6j  /l^  r  f  //  B  ^  ll^  - 


X  ex  i 

and  Tis  the  unique  element  of  X  having  this  property. 

The  minimum  problem  (3/1)  follows  from  equation  (i/2)  when(s)  is 


taken  to  be  the  null  element  of  S. 


REMARK  3. 

The  corresponding  extremal  properties  of  the  interpolating 

—  .  ”7  ^ 

splines  L  Ij^J  are  formally  recovered  by  setting  ^=(?and  B  the 
null  operator. 

S.Z-FINITE  SPACE  Z 

For  the  case  of  specific  interest  here  the  number  of 
"constraints"  is  finite.  The  space  Z  is  consequently  finite 
and,  according  to  thorem  I,  existence  is  automatically 
guaranted  whereas  uniqueness  requires  that  N(T)  be  also  finite. 

Suppose  then  that  ZC^  ,  with  the  usual  inner  product,  and 
that  N(T)  is  of  dimensious  <?  .  Then  L51: 
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C3,/3; 


a)  the  ooerator  A  can  be  expressed  as; 

I)  X  ^  L  /  - - y  ^ 

where  the  \(''are  indeoendent  linear  continuous  functionals  on  X; 

b )  if; 

ir,  J  c?  f? 


t  V  ---  ^ 


then  ; 


fi' t  -  *, 


c)  the  followinq  characterization  theorem  holds  (see  Corollary); 


THEOREM  5.  T  c5  ^  is  the  smoothing  spline  correspcndi  ng  to  (T,  jf, 
Z,  '^  )  iff  there  exist  n  coefficients  U;  such  that; 

I-  K":  ^  N/OTJ 


I 


C~  I 


-  CiC.-yy  ]  ; 


REMARK  ^ 


The  corresoondi nq  theorem  for  interpolating  splines  is  formally 
recovered  by  replacing  the  last  (n)  conditions  in  eqs.  (3, /5)  with 
the  conditions 

Upon  the  remark  I  the  first  equation  characterizes  the  spline 
space  S  corresponding  tc  (T,  Kj).  Hence; 

j3^  J6^=R’" 

and  the  extremal  properties  of  the  smoothing  splines  described  by 
Theorem  4  can  be  formulated  as; 

THEOREM  6 .  If  ^  in  an  arbitrary  element  ofi^andiTis  the  unique 
element  of  S  satisfying  the  condition  (S/fT)  ,  then; 
a )  for  any  x 
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S 


j  A;J  j 


and  any  other  rre  S  having  this  property  belongs  to  the  set: 

-i-hf  LV 

b)  for  any  ':\6S' 

'll 


^  'nu.'yo 

Key 


/ 


y 


(3./6J 


and  r i s  the  unique  element  of  X  having  this  property. 

The  different  classes  of  smoothing  spline  functions  that  can  be 
obtained  from  the  above  abstract  formulation  depend  on  the  choice 
of  X,  Y  'T^and  y 


-38- 


NORMAL  CLOSED  SMOOTHING  SPLINES.  EXISTENCE,  UNIQUENESS  AND  CH ARACTE R I ZA' 


Let  C  be  a  sufficiently  smooth  and  regular  closed  contour.  Take 
and  y-ric)  with  their  standard  inner  products.  Denote 
by  Z  the  curvilinear  coordinate  along  C  measured  from  an  arbitrary 
point  and  normalized  with  resoect  to  the  length  of  C. 

For  the  operator  /  *  ~  tl  take  the  Hi  derivative 

[^with  respect  to  {z.  ).  For  Zc  ^  let  the  functi  onal  s  i/(c‘ ^ 

I 

be  defined  by: 


>  = 


i  1  ^ 


-  L  -- 


with 


^  -  -  —  <C  ^ 

The  null  spaces  of  the  operators  A  andTare  then  given  by /_4 J' 

K'(^J i  ^  ^  I  ^  ;  c-  ei  h'»]  i 

Hence  their  dimensions  are  equal  to  (n)  and  to  (I),  respectvely  so 
that  7/ CTJ/^  orovi  ded  Ti>d  •  The  operator  L  = 

defined  by: 

Z  /  -  x)  A  y J '  Z  Z - y 


so  that : 

The  minimum  problem  (3-1)  reads: 

n.  - 


<  Z'  -n  J 


I  ]^pIz  4^  1 0-44-  ti  7  ' 

'C  ,  r^''  .  k  - 


(^.3) 


x6  yc 

and,  according  to  theorem  I,  it  has  a  unique  solution,  for  any  r, 
as  long  as  vi>  i. 

The  existence  and  uniqueness  ofcK  2r  )  is  thus  established. 

Its  characterization  is  accomplished  through  theorem  3 . 
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F  rom  eq  .  (  3.|5)  with 


o\  and  from  eq.  {4-\)  one  has. 


subsequently • 


ri 

I  £l  ^71^ 


0~ 


•  V>^etn{cJ 


(4.^) 


-i:  -  fL'^,:  -  ; 

where,  according  to  theorem  5,  the^'must  satisfy  the  condition: 


X  e  N  0) 

stating  that  (  1^*0  must  belono  to  )^ij). 

^  As  dimNy/_Tj-i,  eq.  (4-?)  reduces  only  to  the  requirement  that 
^  Kl  be  /^^orthogonal  to  unity.  On  account  of  eq.  {4.5)  this 


C  —  < 

leads  to: 


I 


0 


(3J 


[4.^ 


The  characterization  of  a  function  5*  satisfying  eqs.  (4-^  and(4',i?) 
was  developed  in 

when  dealing  with  closed  interpolating  splines. 

This  finding  is  consistent  with  the  fact  that,  as  mentioned 
in  paragraph  (3),  the  space  of  spline  functions  corresponding  to 
(  T,f4)  is  the  same,  whether  one  deals  with  interpolating  or 
smoothing  splines  [A  I- 

The  details  will  not  be  repeetted  here  and  only  the  relevant 

results  will  be  stated . 

<34.^] 


The  function  ^  can  be  e.xpressed  as: 

:IV_.  -Vl  Q  ,  ^  ^  A: 


i-i  ,  (j-(  jj 

where  the  (n  t^jeoef  f  i  ci  ents  p  ^ '  and  ^Jsati  sfy  the(^  )  equati ons  J 

p-'  r 


JL 
.> 
L- I 

with: 


/  J3i±3i  *  (j -I 


(4  <10) 


(4//) 
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X 

7n  ^  ^ 
L-~  I 


el 


j'  L 

yn  i 


rn. 


r/\  >  o 


yyi  ~  O 


^.‘0  i^j 

EqsVTxpress  the  vanishing  of  the  discontinuities  of  5”  and  of  its 

first  {’^-Z  )  derivatives  at  the  closure  po  i  n  t  cha  racteri  zed  by 

either  values  or  Z-=  I  of^T. 

The  function  given  by  eq.  (4,'^)is  continuous  on  C  to  = 

gether  with  its  first  {^-Z  )  derivatives.  The  discontinuity; 

of  the  )-  tk  derivative  of  c^^'^at  the  point  (9/i  s  given  by; 


The  other  n  equations  needed  to  compute  the  (n+^  )  coefficients 

^  and  ^^are  given  by  eqs.  (<V,6) . 


According  to  eqs,  (i^)  ,  (3,J)  and  (4,^  and  the  definition  of 

innervTh  P  the  space  of  closed  spline^is  defined  by; 


K 


Hence ,  for  any  ><6  hKQ  ^  ^  ^ 

or.T^Iccfunt  Of  eo.  ^  ?£  i; 

This  last  relation  is  the  same  as  that  derived  in  /^i/and  once  again 
reflects  the  fact  that  the  space  Smakes  no  reference  to  the  type 


<yt. 


of  closed  spline,  whether  i nterpol ati ng ^smooth! ng , 

The  closed  smoothing  spline  (T 6  is  obtained  by  integra= 

ting  eq,  q  times  and  using  the  arbitrary  constants  to  impose 

the  continuity  of  S'  and  of  its  first  (  )  derivatives  at  the 

closure  point  so  as  to  make 

The  extremal  properties  formulated  in  section  {-i  )  follow 
from  theoremen  6  with  7*  and  A(:r^^/(,‘5i  ven  by  eq.  (4.^). 


This  concludes  the  proof  of  the  statements  made  in  Section  I  . 
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